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Quark spectra in QCD are linked to fundamental properties of the theory including the identification 
of pions as the Goldstone bosons of spontaneously broken chiral symmetry. The lattice Overlap- 
Dirac operator provides a nonperturbative, ultraviolet-regularized description of quarks with the 
correct chiral symmetry. Properties of the spectrum of this operator and their relation to random 
matrix theory are studied here. In particular, the predictions from chiral random matrix theory in 
topologically non-trivial gauge field sectors are tested for the first time. 



PACS numbers: 11.15.Ha, 11.30Rd, 12.38Gc 

An important property of massless QCD is the sponta- 
neous breaking of chiral symmetry. The associated Gold- 
stone pions dominate the low-energy, finite-volume scal- 
ing behavior of the Dirac operator spectrum in the micro- 
scopic regime, defined by 1/Aqcd << L « I/ttItt, with 
L the linear extent of the system ^. This behavior, in 
turn, can be characterized by chiral random matrix the- 
ory (RMT), which lead to a revival of RMT, first used to 
understand the energy levels of nuclear matter [g| . What 
enters into the RMT description of the low-energy, finite- 
volume scaling behavior are some symmetry properties 
of the Dirac operator and the sector of fixed topological 
charge under consideration [^l). The RMT predictions 
are universal in the sense that only the symmetry prop- 
erties, but not the form of the potential matters [g. 
Furthermore, the properties tested in this letter can be 
derived directly from the effective, finite-volume parti- 
tion functions of QCD of Leutwyler and Smilga, without 
the detour through RMT §, though RMT nicely and 
succinctly describes and classifies all these properties. 
The topological charge enters the RMT prediction via 
the number of fermionic zero modes, related to the topo- 
logical charge through the index theorem. The symme- 
try properties of the Dirac operator fall into three classes, 
corresponding to the chiral orthogonal, unitary, and sym- 
plectic ensembles B . Examples are fermions in the fun- 
damental representation of gauge group SU(2) for the 
orthogonal ensemble, fermions in the fundamental rep- 
resentation of gauge group SU(-/Vc) with Nc>i for the 
unitary ensemble, and fermions in the adjoint representa- 
tion of gauge group SU(-/Vc) for the symplectic ensemble. 

The classification according to the three different chi- 
ral random matrix ensembles is intimately connected to 
the chiral symmetry properties of the fermions. A good 
non-perturbative regularization of QCD should therefore 
retain those chiral symmetry properties. Until recently 
such a regularization was not known. The next best 
thing were staggered fermions, which at least retained 



a reduced chiral-like symmetry on the lattice. Indeed, 
staggered fermions were used to verify predictions of chi- 
ral RMT R,H|, albeit with two important shortcomings: 
(i) staggered fermions in the fundamental representation 
of SU(2) have the symmetry properties of the symplectic 
ensemble g], not the orthogonal ensemble as continuum 
fermions, while adjoint staggered fermions belong to the 
orthogonal ensemble [Q, not the symplectic one. (ii) 
staggered fermions do not have exact zero modes at fi- 
nite lattice spacing [Q , even for topologically non-trivial 
gauge field backgrounds, and thus seem to probe only the 
V ~ Q predictions of chiral random matrix theory |0,|[ . 

The development of the overlap formalism for chiral 
fermions on the lattice ||l2[ recently lead to the mass- 
less Overlap-Dirac operator, a lattice regularization for 
vector-like gauge theories that retains the chiral proper- 
ties of continuum fermions on the lattice p3| . In par- 
ticular, the continuum predictions of chiral random ma- 
trix theory should apply. Overlap fermions have exact 
zero modes in topologically non-trivial gauge field back- 
grounds |14|, allowing, for the first time, verification of 
the RMT predictions va v ^ sectors. The nice agree- 
ment we shall describe further validates the chiral RMT 
predictions on the one hand and strengthens the case for 
the usefulness of the Overlap regularization of massless 
fermions on the other hand. 

The massless Overlap-Dirac operator [pSi is given by 



D=]^[l + -i^e{H^{m))] 



(1) 



Here, 75iJ^(— to) is the usual Wilson-Dirac operator on 
the lattice and e denotes the sign function. The mass m 
has to be chosen to be positive and well above the critical 
mass for Wilson fermions but below the mass where the 
doublers become light on the lattice. In this letter, we 
will be interested in the low lying eigenvalues oi H = 
75-D, which is a hermitian operator. Relevant properties 
of this operator can be found in Ref . (lil . We will use the 
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FIG. 1. Plots of the unfolded level spacing distribution P{s) versus s for the various ensembles. The three curves are 
the predictions for the three ensembles where the sharpest curve is for the symplectic ensemble, followed by the curve for the 
unitary ensemble and the broadest curve is for the orthogonal ensemble. 



Ritz algorithm ||l^ applied to H^ to obtain the lowest few 
eigenvalues. The numerical algorithm involves the action 
of iJ on a vector and for this purpose one will have to use 
a representation of e{H^{m)). We will use the rational 
approximation discussed in Ref. ||l4| , [l6|] . 

For fermions in the fundamental representation of 
gauge group SU(2), the Overlap-Dirac operator is real 
since the underlying Wilson-Dirac operator is real [[l7[ . 
They are therefore expected to fall into the chiral orthog- 
onal ensemble. For fermions in the adjoint representation 
the spectrum of the Overlap-Dirac operator is doubly de- 
generate, since the spectrum of the underlying Wilson- 
Dirac operator is doubly degenerate |l8[, and they are 
therefore expected to belong to the symplectic ensemble. 
For three or more colors the Overlap-Dirac operator in 
the fundamental representation is complex with no addi- 
tional symmetries and is therefore expected to fall into 
the chiral unitary ensemble.n 

We will first present our results for the distribution of 
the "unfolded" level spacing Mv^ ■ The unfolding of the 
eigenvalues of H are done in the following manner. Let 
E" label the non-zero positive eigenvalues of H with n 
labeling the configuration number and £^" > E^_i for 
all i. We are considering only the positive eigenvalues 
of H, since the non-zero eigenvalues of H all come in 
positive/negative pairs [ |l4| . For fermions in the adjoint 
representation of SU(iVc) the spectrum of H is doubly de- 
generate and we only keep half the spectrum, dropping 
the degeneracy. Unfolding proceeds by first sorting all 
the E" in ascending order and associating the location 
iVf of El" in the sorted hst with El". N^ is referred to 
as the unfolded spectrum and the level spacing is sim- 
ply given by (A^,'+i — N")/N where N is the number of 



* Fermions in U(l) abelian gauge theory are another example 
of the unitary case. For numerical tests in two dimension, see 



configurations. The distributions of the unfolded level 
spacing, s, in RMT are well approximated by the various 
Wigner distributions M 



Pi.s) = 



3s\- 



262144 „4„-#^s 



64 2 



7297r' 



xs e 



orthogonal ensemble 
unitary ensemble 
symplectic ensemble 



(2) 



In our numerical simulations we computed the low ly- 
ing spectrum of the Overlap-Dirac operator in the fun- 
damental representation on pure gauge SU(2) configura- 
tions with (3 — 1.8 as an example of the chiral orthog- 
onal ensemble, on pure gauge SU(3) configurations with 
/3 = 5.1 as an example of the chiral unitary ensemble, and 
in the adjoint representation on pure gauge SU(2) con- 
figurations with (3 = 2.0. The lattice size was 4^ in all 
cases. The various level spacing distributions are shown 
in Fig. ^ There is very clear evidence that the SU(2) and 
SU(3) ensembles with fermions in the fundamental repre- 
sentation fall into the orthogonal and unitary ensemble, 
respectively, and the SU(2) ensemble with fermions in the 
adjoint representation falls into the symplectic ensemble. 

We next turn to the distribution of the lowest eigen- 
value for the various ensembles. Chiral RMT predicts 
that these distributions are universal when they are clas- 
sified according to the three ensembles and according to 
the number of exact zero modes v within each ensemble 
and then considered as functions of the rescaled variable 
z — SFAmin- Here V is the volume and S is the infinite 
volume value of the chiral condensate {'ipip) determined 
up to an overall wave function normalization, which is 
dependent in part on the Wilson-Dirac mass m. RMT 
gives for the distribution of the rescaled lowest eigen- 
value for the orthogonal ensemble, expected to apply to 
the fermions in the fundamental representation of SU(2), 
in the v — and i^ = 1 sector |2l| 
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FIG. 2. Plots of Pmin{z) vei'sus 2 for the various ensembles in the lowest two topological sectors. The curve in each plot is 
a fit to the prediction from random matrix theory with the best value for the chiral condensate. 



For the unitary ensemble, expected to apply to the 
fermions in the fundamental representation of SU(iVc) 
with Nc > 3, the RMT predictions are [|l|J22| 



-* mill V^J 



(4) 

Finally, for the symplectic ensemble, expected to apply 
to the fermions in the adjoint representation, the RMT 
prediction is il|,|3|,§ 
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We are interested in i^ = 1 since the eigenvalues are 
doubly degenerate. A closed form expression is not 
known for T^{z'^). Instead, a rapidly converging se- 
ries [^,|1 based on partitions of integers is available, 
namely T^{x) = 1 + X^dLi o.d.x'^ where 
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{2v + 2j - i + 4) ' 
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(6) 



Here, the integer partition k — {ki,K2, . . . ,Kd} has 
length d and weight 1(k,) — "^^Ki- A pair of integers 
is associated with n — {{i,j)\l < i < ^k),1 < j < k.^}, 
and Kf = CaTd{j\Kj > i) is the conjugate partition. 

We compare the RMT predictions with our data in 
Fig. S. If one knows the value of the chiral condensate 
in the infinite volume limit, E, the RMT predictions for 
Pniin(z) are parameter free. On the rather small systems 
that we considered here, we did not obtain direct esti- 
mates of S. Instead, we made one-parameter fits of the 
measured distributions, obtained from histograms with 
jackknife errors, to the RMT predictions, with S the free 
parameter. Our results, together with some additional 
information about the ensembles, are given in Table |. 
We note the consistency of the values for E obtained in 
the z^ = and 1^ = 1 sectors of each ensemble. Alter- 
natively, we could have used the value of E obtained in 
the 1/ = sector, to obtain a parameter free prediction 
for the distribution of the rescaled lowest eigenvalue in 
the 1/ = 1 sector. Obviously, the predictions would have 
come out very well. 

For the two ensembles with the fermions in the fun- 
damental representation, we also found 81 (for SU(2)), 
and 147 (for SU(3)) configurations with two zero modes 
and 1 and 3 with three zero modes. For the orthogonal 



TABLE I. The chiral condensate, E, from fits of the distri- 
bution of the lowest eigenvalue to the RMT predictions. The 
third column gives the Wilson-Dirac mass parameter used, 
the fourth the number of configurations, A^^ , in each topolog- 
ical sector and the last the confidence level of the fit. 
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254 


0.2900(30) 
0.2931(45) 
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0.313 



ensemble, we are not aware of a prediction for Pmin(^) in 
the V = 2 sector, while for the unitary ensemble our data, 
albeit with very limited statistics, agrees reasonable well 
with the parameter free prediction, eq. (jj) with S from 
Table |. 

For fermions in the adjoint representation, we keep 
only one of each pair of doubly degenerate eigenvalues, so 
I' = 1 is the sector where there are two exact zero modes. 
Such gauge field configurations cannot be assigned an in- 
teger topological charge since integer charges give rise to 
zero modes in multiples of four p4], and we note there 
are a significant number of configurations with two zero 
modes as seen in Table H. The good agreement with the 
RMT prediction found in this case lends further support 
to the existence of configurations with fractional topo- 
logical charge |Q. 

In this letter we have tested the predictions of chiral 
random matrix theory using the Overlap-Dirac operator 
on pure lattice gauge field ensembles. We find full agree- 
ment with the unfolded level spacing distributions on all 
three ensembles. We also found that the distribution of 
the lowest eigenvalue in the different topological sectors 
fitted well with the predictions of chiral RMT, with com- 
patible values for the chiral condensate from the different 
topological sectors. This is the first test of the influence 
of topology on the Dirac spectrum in the microscopic 
regime. 
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